A criterion for entanglement in two two-level systems 
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We prove a necessary and sufficient condition for the occurrence of entanglement in two two-level 
systems, simple enough to be of experimental interest. Our results are illustrated in the context of 
a spin star system analyzing the exact entanglement evolution of the central couple of spins. 
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It is today experimentally possible to realize many en- 
tangled states in different contexts like CQED , 
metal-superconductor junctions ||, @ and other solid 
state systems Q, @|. A very debated and open issue 
concerns the individuation of efficient methods to detect 
and quantify entanglement. In principle, tomography 
techniques allow the reconstruction of the state of the 
system and to establish whether the system is entangled 
with the help of some separability criteria. To deter- 
mine in laboratory the density operator matrix elements 
requires however considerable efforts, making thus desir- 
able to probe and to quantify the occurrence of entan- 
glement without using a tomographic approach. To this 
end Bell inequalities and witness operators [§)-[l2j could 
be used even if their exploitation is indeed still an hard 
task. It is indeed an exciting challenge to find criteria 
for detecting entanglement involving the measurement 
of few observables of clear physical meaning. Methods 
based on violation of local uncertainty relations [l3| or 
inequalities for variances of observables [lij have been 
quite recently proposed. In this paper we present a crite- 
rion for entanglement well suited for bipartite systems in 
which each component is a two-level system. Our results 
do provide an incisive tool to infer the existence of non 
classical correlations from few and simple measurements. 
Our criterion is successfully applied to a spin star system. 

Consider a bipartite system composed by two two- 
level systems described by the Pauli operators <jj = 
(o± , a z 1 ^ ) i = 1,2 respectively. Suppose that its den- 
sity operator, in the factorized basis {|Tt), |T4), lit); 
| 11}} of the eigenstates of has the following 

structure 
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The quite simple form of /?, given by eq.([T]) naturally 
arises in many physical scenarios not necessarily involv- 
ing spin ~ systems @, d, [H, [Hj]. Since the concurrence 
function associated to eq.([T]) is 



Cone = max[0, 2(|c| — \/ae)] 



(2) 



the existence of entanglement in our system requires that 
the geometric mean of the two populations a and e be- 
comes less than the coherence amplitude \c\. The two 



probabilities a and e of finding the system in the state 
t|) and | J, J,) respectively, determine the mean value of 



• <Jz ) and as follows 
(S z ) = tx{pS z } = a- 



(3) 



((S z ) 2 )^t r {p(S z ) 2 } = a + e (4) 
which in turn imply that 

(AS z f = ((S z f) - ((S z )f = (a + e) - (a - e) 2 . (5) 
In accordance with the Landau's condition 11611. 



\c\ < Vbd, 



(6) 



we may thus state that the presence of entanglement in 
the two two- level systems described by eq.([T]), necessarily 
requires 



me < Id < v bd. 



(7) 



Exploiting such an inequality in eq.^ yields the follow- 
ing upper limit to the variance of S z : 



(AS Z ) 2 < Abd- {b + d)[(b + d) - 1]. 



(8) 



It is of relevance the fact that when the density matrix 
(U) assumes the special form 



Pit) 



(a 

b b 

b b 

\ e 



(9) 



the necessary condition JS} simply reduces to 



or, equivalcntly, 



(AS Z ) 2 < 2b 



((S z )) 2 > 1 - 46. 



(10) 



(11) 



Since when condition (fTTj) is fulfilled then the concur- 
rence is different from zero, we may claim that eq. (jlip 
is a necessary and sufficient condition in order to have 
entanglement in two two-level systems described by p 



2 



as given by eq.©. Stated another way the occurrence 
of entanglement in the system may be checked simply 
comparing the square of the mean value of S z with the 
population b. It is however important to stress that the 
result we have obtained allows us to immediately deduce 
that if b > i , then we can be sure that the two two-level 
systems are entangled. On the contrary when b < j we 
need also to know the mean value of S z . In this case 
indeed the bipartite system is entangled if, and only if, 

- \2b - 1| < (S z ) < -VT^4b (12) 

or, alternatively, 

VT T7 46< (S z ) < \2b- 1|. (13) 

In other words when b < i, the presence of entanglement 
in the system is compatible only with S z mean values 
'squashed' toward the extremes of its interval of variabil- 
ity [—1,1]. In turn it implies that under the condition 
< b < \, in order to have entanglement, the probability 
of finding the system in the state ||T) must be apprecia- 
bly different from the probability to find the system in 
the state \[[). 

Our criterion may be successfully exploited in order 
to analyze the entanglement evolution in a spin system 
describing a physical scenario of interest in many physical 
contexts. Consider indeed a system constituted by two 
uncoupled spins ^, denoted by A and B and hereafter 
called central system, each one interacting with M — 2 
mutually uncoupled spins ^. In particular we suppose 
that the two central spins interact with each of the M — 2 
spins in the same way. This system, known as spin star 
like system [l?], OH, can be described by adopting the 
following hamiltonian model 

H = Hq + Hj (14) 

with 

H =uj(S z + J z ), (15) 



Hi = a{S+J- +S-J+), (16) 



where S z and S± = \{o^± + c±^) are spin operators 
acting in the Hilbert space of the central system and J z 
and J± are the collettive spin operators describing the 
other M — 2 spins. This hamiltonian model can be suc- 
cessfully used to describe for example electronic spins 
in semiconductor quantum dot coupled by hyperfine in- 
teraction with nuclear spins, or electronic spins bound 
to phosphorus atoms in a matrix of silica or germanio in 
presence of defects [lj| . The hamiltonian model (fT4")l pos- 
sesses permutational symmetries successfully exploitable 
to exactly solve the relative time-dependent Schrodinger 
equation [2(|, [2l[. Suppose that at t = the central 
system is in a common eigenstate of S* 2 = (Si + 52 ) 2 
and S z denoted by IS*, Ms). At the same time the re- 
maining M — 2 spins are supposed in the state | J, Mj, v) , 
eigenstate of the collettive angular momentum operators 
J 2 and J z . The index v, depending on J, allows us to 
distinguish between different states of the coupled angu- 
lar momentum basis characterized by the same J and 
Mj. Hamiltonian (|14[) is invariant by permutation of 
the two central spins as well as of an arbitrary couple of 
spins among the M — 2 of the second block. Moreover 
[S 2 ,H] = [J 2 ,H] = [S z + J Z ,H} = [Jf nt ,H] = 0, J mt 
being an intermediate angular momentum resulting from 
the coupling of selected at will individual angular mo- 
mentum of the M — 2 spins. At a generic time instant t 
we can write 



\rf,(t)) = e- lHot e- lH ' t \S,M s )\J,Mj,v) 1 (17) 



being [H ,Hi] = 0. The case S = and consequently 
Mg = 0, is a trivial one corresponding to an eigenstate 
of the Hamiltonian whatever the value of J and Mj are. 
When instead S — 1 we have 



Hf n \l,M s )\J,Mj,v) 



+ r Ms \l,-M s )\J,Mj + 2M S , v) ) 



p Ms \1,M s )\J,Mj,v)+ 



(18) 



and 



Hf n+L \l,M s )\J,Mj,y) = 

_ j a 2n +HV2) 2n+1 PMs(P 2 M s +r 2 Ms ) n \l,0)\J,Mj + M s ,iy) 



if M s ^0, (19) 



2 ™ +1 (V2) 2 " +1 [p-^p^+plT^W.Mj-l^) +p 1 (p 2 _ 1 +p 2 )"|l,-l)|J,M J + l, l /)] if Ms -0 



with 



Ps = y/J(J + l)-Mj(Mj + s), s = ±1 

Po = 1, 



r s = VJ(J + 1) - Wj + S )( M J + 2s ). a = ±1 
r = 0. 



(21) 



(20) 
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After straightforward calculations we thus get 

|V>(t)) = e --(M s +M j)t ^ Mj {t)\l,M s )\J,Mj,v 
+ Bf^(t)\l, -M s ) | J, Mj + 2Af s , i/) 

- *<W o [^flr, (t) 1 1, 1) | J, Mj - 1, u)+ 
+E^ j (t)\l,-l)\J,M J + l,p)]}. 



with 
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r 2 
r M s 

~r 2 I 2 

PjVf s + r M s 



(22) 



3 (v /2 ^M S + r M S )^) + 



if Ms = — 1. Suppose to prepare the two spins of the 
central system and the remaining M — 2 = N spins in 
the state 



|V(0)) = \l,l)\N/2,-N/2 + k,\) 



(31) 



(23) 



where k = 0, 1, — N. Thanks to our results expressed by 
cq. (flTJ)) . we can state with certainty that at a generic time 
instant f the two central spins are entangled if and only 
if (AS z f < 26(f). In Fig.Q we plot (AS Z ) 2 and 26(f) 
in correspondence to N = 100 and k = 2 as a function of 
at. Since condition (|10p is not verified whatever the time 
instant f is, we can easily conclude that starting from the 
initial condition (|3ip with N — 100 and k = 2, the central 
system is unable to develop quantum correlations. At the 
light of the result before discussed, the physical reason 
of this incapacity of the system to develop entanglement 
stems from the fact that the time evolution from this 
initial condition never makes significantly different the 
probabilities of finding the two spins in the state Iff) or 
Iff) respectively. 



pMs f.\ _ PMs r M s 
n JM.,\ l ) - 2 T r 2 

VMs ' Ms 



Df Mj (t) = 



PMs 



PMs ^ T M S 



(cos(^2(p 2 Us +r 2 Ms )at) - l), 

(24) 



: sin 



J M S ^ M, 



P-i 



P-i +Pi 



Pi 



P-i +P 2 i 



(25) 

S m^2(pl 1+ pl)aty (26) 
: sin( v /2(^ + P \)at). (27) 



We are interested in the entanglement dynamics of the 
two central spins A and B. Tracing p(t) = |i/ , (f))(^(f)|, 
over all the degrees of freedom of the M — 2 spins around 
A and B we get the reduced density matrix of the central 
system. It assumes the block diagonal structure given by 
eq.@, every matrix element being an explicit function 
of the probability amplitudes (|23p - ([?71) . In particular if 
M s = 1 

a(t) = K| fj (f)| 2 , 6(f) = |C^ f/ (f)| 2 , e(f) = K| f/ (f)| 2 , 

' (28) 

while if M s = 

*(*) = IdYmM 2 ' W = IAjmM 2 + Km. 7 WI 2 , 

e(f) = |^(f)| 2 



(29) 



and 



<*(*) = Km.WI 2 , Ht) = |Cf^(f)| 2 , e(f) = K| f; (f)| 2 

' (30) 
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FIG. 1: {AS Z ) 2 (dashed line) and 26(f) (bold line) against 
at in correspondence to k = 2 and N — 100 




FIG. 2: (AS Z ) 2 (dashed line) and 26(f) (bold line) against at 
in correspondence to k = 98 and N = 100 

It is of relevance the fact that the same conclusion is 
reached starting from |^(0)) = |1, l)\N/2,-N/2 + k, 1) 
with k < 50. The behaviour of the system is instead 
deeply different when k exceeds 50. As an example we 
compare in Fig.© {AS Z ) 2 and 26(f) for k = 98 and 
N = 100. In this case, there exist different time intervals, 
in which the two spins A and B are entangled because 
(AS^) 2 < 26(f). Thus the parameter k controlls the abil- 
ity of the system to generate entanglement in the central 
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system. We expect that the amount of entanglement also 
depends on the choice of k that in turn determines the S z 
fluctuations. This indeed is true as witnessed by Figures 
© and (|1]) where the time evolution of the concurrence 
is reported. We stress that for k = 100 the initial con- 
dition is an exact stationary state while for k = 99 the 
two central spins reach the maximum compatible degree 
of entanglement since in this case the population of \H) 
exactly vanishes at any time instant so that eq. tfTOf is 
always fulfilled except when b(t) = 0. On this basis we 
foresee and we have proved a specular behaviour starting 
from the initial condition |1, -l)\N/2, -N/2+k, 1) in the 
sense that in this case the entanglement reaches its max- 
imum when k — 1. We emphasize that the constraint 




FIG. 3: Concurrence function of the central system 
against at and k in correspondence to the initial condition 
1, l)|iV/2, -N/2 + k) for N = 100 




FIG. 4: Concurrence function of the central system against 
at and k (97 < k < 100) in correspondence to the initial con- 
dition [1, 1) |JV/2, -N/2 + k) for N = 100 

on the fluctuations of S z expressed by ea. (fT0)) is the key 
of the simplicity with which we analyze the appearance 
and the disappearance of entanglement as a function of 
time and more important to understand its dependence 
on the value of k. In this paper we have proved that 
when two two-level systems are described by a density 
matrix expressed by eq.((9]) at any time instant t, there 
is a simple and reliable procedure of experimental inter- 
est after which the occurrence of quantum correlations 
be surely claimed or excluded. We propose indeed the 
measurement of at most two populations which amounts 
at comparing the variance of S z with the probability of 
finding one spin up and the other down. Our approach 
to the entanglement of the pair of two level systems has 
the merit of directly involving quantities having a clear 
physical meaning. Applying our criterion to a spin star 
system we are able to fully exploit the novelty of our 
point of view to explain the dependence of the ability of 
the system to develop entanglement on the initial condi- 
tions. 
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